Soliton trains in Bose-Fermi mixtures 
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We theoretically consider the formation of bright soiitons in a mixture of Bose and Fermi degen- 
erate gases. While we assume the forces between atoms in a pure Bose component to be effectively 
repulsive, their character can be changed from repulsive to attractive in the presence of fermions 
provided the Bose and Fermi gases attract each other strongly enough. In such a regime the Bose 
component becomes a gas of effectively attractive atoms. Hence, generating bright soiitons in the 
bosonic gas is possible. Indeed, after a sudden increase of the strength of attraction between bosons 
and fermions (realized by using a Feshbach resonance technique or by firm radial squeezing of both 
samples) soliton trains appear in the Bose-Fermi mixture. 
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Solitonic solutions are a very general feature of nonlin- 
ear wave equations. Soiitons have been studied in many 
different physical systems ranging from particle physics 
to optics. They differ from ordinary wave packets as they 
retain their shape while propagating instead of spreading 
due to dispersion. This intriguing feature is based on the 
existence of a nonlinear interaction which compensates 
for dispersion and produces a self-focusing effect on the 
propagating wave packet. 

Dilute atomic quantum gases offer a unique environ- 
ment to study fundamental solitonic excitations in a pure 
quantum system with intrinsic nonlinearity. Since the 
interparticle interaction causing this nonlinearity can be 
both attractive and repulsive, the Gross-Pitaevskii equa- 
tion describing the evolution of the condensate wave func- 
tion exhibits both dark and bright solitonic solutions pj . 
Dark soiitons as a fundamental excitation in stable Bose- 
Einstein condensates with repulsive interparticle interac- 
tion have been studied in different geometries 
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Bright soiitons have been observed in Bose-Einstein 
condensates of 7 Li in quasi-one-dimensional geometry 
[HEi- However, in three-dimensional geometry usually 
used to prepare the sample the necessary large and neg- 
ative scattering length leads to density-limited particle 
numbers (dynamical instability - collapse) . The observa- 
tion of bright soiitons was therefore only possible due to 
magnetic tuning of the interactions from repulsive (used 
to form a stable Bose-Einstein condensate) to attractive 
during the experiments. 

Another experimental approach to bright matter wave 
soiitons was realized in the recently reported observa- 
tion of gap soiitons in a condensate with repulsive in- 
teractions by engineering of the matter wave dispersion 
relation via sophisticated manipulation in aperiodic po- 
tential (concept of negative effective mass 

In this Letter we propose a novel scheme to realize 
bright soiitons in one-dimensional atomic quantum gases. 
In particular we study the formation of bright soiitons in 
a Bose-Einstein condensate embedded in a quantum de- 
generate Fermi gas. One important feature is that this 
mixture allows tuning of the one-dimensional interac- 



tions not only by Feshbach resonances but also by simply 
changing the trap geometry. 

We consider the bare interaction between bosonic 
atoms to be repulsive (gs > 0) whereas the particle 
interaction between bosonic and fermionic atoms is as- 
sumed to be strongly attractive (gBF < 0). Bright soli- 
tons in Bose-Fermi gas mixtures are then produced as 
a result of a competition between two interparticle in- 
teractions: boson-boson repulsion versus boson-fermion 
attraction. Experimentally, this situation is for example 
accessible in Bose-Fermi mixtures of bosonic 87 Rb atoms 
and fermionic 40 K atoms. 

We determine the critical strength of attraction be- 
tween bosons and fermions necessary for the forma- 
tion of bright soiitons and show that these parameter 
regimes might be achievable in present experiments. We 
study the formation of bright soiitons by different ex- 
citation mechanisms: either increasing the attractive 
boson-fermion interaction by Feshbach resonance tech- 
niques or by radial squeezing of the mixture which corre- 
sponds to increasing the effective one-dimensional scat- 
tering length. We contrast the response of the system fol- 
lowing adiabatic and fast increase of the boson-fermion 
interaction strength. 

We consider a Bose-Fermi mixture confined in a trap 
at zero temperature and describe this system in terms of 
the many-body wave function ^(xi, xn b ; yi, yjv F ), 
where Nb and Nf are the numbers of bosons and 
fermions, respectively. However, instead of direct 
(and of course approximate) solving of the many-body 
Schrodinger equation we start with the equivalent ap- 
proach based on the Lagrangian density. Since the 
fermionic sample is spin-polarized only boson-boson and 
boson-fermion interactions are included. At zero tem- 
perature we assume that the wave function of Bose-Fermi 
mixture is a product of the Hartree ansatz for bosons and 
the Slater determinant (antisymmetric wave function) for 
fermions 
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teract. Then we have 
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Now, all the single-particle wave functions (ip^ B \ 

ip\ ', -•■,<Ppf ] l) have to be determined. Therefore the 
many-body wave function Q is inserted into the La- 
grangian density and integrated over the spatial coor- 
dinates to get the Lagrangian. In fact, this integration is 
performed only over Nb — 1 bosonic and N F — 1 fermionic 
coordinates and leads in this way to the mean-field single- 
particle Lagrangian. The Euler-Lagrange equations cor- 
responding to this Lagrangian are the basic equations of 
the presented approach 
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All of the above equations have simple interpretation. 
Removing the last terms in these equations (i.e., neglect- 
ing the mean-field interaction energy between Bose and 
Fermi components in comparison with other energies) 
one recovers the Gross-Pitaevskii equation for degener- 
ate Bose gas and the set of Schrodinger equations de- 
scribing the noninteracting Fermi system. It is easy to 
notice that when the Bose and Fermi gases attract each 
other strongly enough the mean-field energy connected 
with this attraction can overcome the repulsive mean- 
field energy for bosons. This can happen provided both 
phases have a significant overlap. It means that the pres- 
ence of a degenerate Fermi gas changes the character of 
the interaction between bosonic atoms from repulsive to 
attractive. Therefore the formation of bright solitons be- 
comes possible. 

Based on the above considerations one can easily write 
the condition for the value of critical strength of attrac- 
tion between bosons and fermions. It is given by gsn-B = 
\gg F \n F , where n B = N B W [B) \ 2 and n F = |^| 2 
are the densities of both fractions, normalized to the 
number of particles, taken at the center of the trap. 
A rough estimation of g c g F can be found assuming the 
densities of both components are calculated within the 
Thomas- Fermi approximation and components do not in- 
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where C = d (a B /a B ) 3/5 K/af ) 3 X^/X]/ 2 , d = 
39/10 52/5 1 a± [ s tj ie radial harmonic oscillator 

length, A = uj z /uj± defines the aspect ratio of axially sym- 
metric trap, and as is the s-wave scattering length for the 
pure Bose gas related to the interaction strength through 
cjb = 47rft 2 as/TOB. The condition © has several impli- 
cations. Squeezing radially both Bose and Fermi compo- 
nents decreases the value of critical gsF- The same hap- 
pens when the number of fermions is getting bigger in 
comparison with the number of bosons. For a particular 
trap numbers of atoms are limited by the occurrence of 
a collapse [j| . In the case of experiment of Ref. jf} it was 
found that the system was stable if the number of atoms 
in both species ( 87 Rb and 40 K in |2, 2 > and |9/2, 9/2 > 
hyperfine states, respectively) were smaller than 2 x 10 4 . 
Taking parameters of that experiment (uj b = 2ir x 215 
Hz and uif = 2ir x 16.3 Hz) and assuming the following 
numbers of atoms Nb = 10 3 and N F — 10 4 , one gets the 
critical coupling \gg F \ = 6.7 gs which equals the natu- 
ral value of gsF for 87 Rb - 40 K mixture in the double- 
polarized state mentioned above. It is understood that 
the relation 10 is only the necessary condition for cre- 
ation of bright solitons. Another important factor is the 
geometry of the system. 
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FIG. 1: Density profiles (normalized to one) of one- 
dimensional Bose-Fermi mixture consisting of Nb = 1000 
bosons and N F — 100 fermions for various strengths of at- 
traction between components as described in the label. The 
effective repulsion for bosons is governed by gs = 0.0163. 
Solid lines correspond to the Bose fraction whereas the dashed 
ones indicate fermions. By increasing the attraction between 
components some fermions are pulled inside the Bose cloud. 
All quantities are given in oscillatory units calculated based 
on fermionic component. 

In the following we will concentrate on one-dimensional 
geometry, which is most favorable for the appearance 
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of solitons. In Fig. ^ we plot the density profiles of 
Fermi (dashed lines) and Bose (solid lines) components of 
one-dimensional mixture in its ground state. The gases 
are confined in traps with frequencies 2n x 30 Hz (for 
fcrmions) and 2tt x 20 Hz (for bosons). To get these 
curves we solve numerically the set of Eqs. by evolv- 
ing adiabatically the coupling constant gsF from zero to 
the given value. Even though the number of bosons is 
10 times larger than the number of fermions, due to the 
Pauli exclusion principle the size of the fermionic cloud is 
much bigger. We see that after turning on the attractive 
forces between components some number of fcrmions is 
drawn inside the bosonic cloud. When the attraction is 
increasing further the fermionic cloud is clearly divided 
into two distinguishable parts. One of them is a broad 
background gas whereas the second is the narrow den- 
sity peak hidden within the bosonic peak. Both peaks 
get higher and narrower when the attraction is getting 
stronger. This is a sign of effective attraction. 
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FIG. 2: Illustration of the solitonic character of the ground 
state of Bose-Fermi mixture. Both frames show the densities 
34 ms after switching off the trapping potential. Only when 
the attraction between the species is stronger than the critical 
one (the lower frame) the system forms double-peak structure 
which does not spread in time. 

It turns out that for strong enough attraction between 
Bose and Fermi components the central peaks (bosonic 
and fermionic ones) in Fig. ^form, in fact, the soliton. 
After switching off the trapping potential the fermionic 
background is lost but the bosonic peak persists without 
changing its shape and confines fermionic density whose 
profile is also preserved. It is illustrated in Fig. [3 where 
fermionic and bosonic densities are plotted some time af- 
ter the opening the trap. Upper and lower frames differ 
by the value of gBF ■ In the case of the upper frame 
the strength of attraction is weaker than the critical 
one and no solitonic behavior is observed - the bosonic 



and fermionic clouds spread out. For the lower frame 
the attraction is strong enough and the formation of 
a bright soliton is observed. Such a structure can be 
forced to move by imposing a momentum on it (realized 
experimentally by applying the Bragg scattering tech- 
nique) and its shape again does not change. The critical 
value of the coupling constant is in our one-dimensional 
model approximately equal to g^ F ~ —1.0 osc. units 
and can be compared to the value obtained based on 
the one-dimensional counterpart of condition (0. Tak- 
ing from Fig. ^the ratio nB(0)/np(0) = 57.5 one gets 
9b f = — 0.94 osc. units what remains in agreement with 
numerical estimation. 
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FIG. 3: Density profiles of one-dimensional Bose-Fermi mix- 
ture after switching the attraction between Bose and Fermi 
components from gsF = —0.5 to qbf = —2.0 oscillatory 
units. The interaction strength is changed linearly during 
8.5 ms (the upper frame) and 17 ms (the lower frame). The 
snapshots are taken at 20 ms. More adiabatic change of the 
strength of attraction results in lower number of solitons. 

In Figs. |3 and 0] we show the densities of Bose and 
Fermi components after the strength of attraction be- 
tween bosons and fermions has been increased by us- 
ing the Feshbach resonance technique. The other way of 
changing the interaction strength could be the firm ra- 
dial squeezing of both samples. The basic observation is 
that after some time after switching the mutual interac- 
tion the bosonic cloud breaks into several peaks provided 
the final coupling (attraction) between both components 
is strong enough. Stronger attraction results in bigger 
number of peaks. Each bosonic peak (marked by the 
solid line) contains the fermionic density (marked by the 
dashed line). Such double-peak structures oscillate al- 
most without changing their shapes. Switching off the 
trap shows that the observed structures are indeed soli- 
tons. In the case of Fig. 0]the strength of the coupling be- 
tween Bose and Fermi gases is changed instantaneously. 
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The contrast of the peaks depends on the initial width 
of the bosonic cloud and it gets higher for bigger number 
of bosons or smaller bosonic trap frequency. The latter 
can be realized only in the optical trap. On the other 
hand, in Fig. [3]we show the response of the system when 
the mutual attraction is increased within the finite time 
which is of the order of trap period. Slower switching of 
the interaction leads to a smaller number of solitons. 
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FIG. 4: Density profiles of one-dimensional Bose- Fermi mix- 
ture 10.5 ms after changing the strength of attraction between 
fractions from the initial value qbf = —0.8 osc. units to the 
final one indicated by the label. Here, the strength of the 
interaction is increased instantaneously. The bosonic cloud 
breaks into several bright solitons, each soliton contains a 
piece of fermionic cloud. The bigger jump in the coupling 
<jB f results in a bigger number of solitons. 

Based on numerical results, we propose two schemes 
for generating bright solitons in degenerate Bose-Fermi 
mixtures. First of all, the system has to be pushed in 
the range of strong enough attraction between fermions 
and bosons. This can be achieved by producing a mix- 
ture in an appropriate doubly-polarized state (for details 
see E3) and by using the Feshbach resonance technique. 
Another way is to follow the idea reported in Ref. [Tij ac- 
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cording to which squeezing radially both samples firmly 
enough imposes a one-dimensional geometry on the sys- 
tem with one-dimensional interaction strength given by 
<73-D/ 7ra i (a± is the radial harmonic oscillator length). 
Increasing radial confinement leads then to higher effec- 
tive one-dimensional interaction strength. Therefore, one 
way to generate solitons in Bose-Fermi mixture would fol- 
low the scenario: (a) the ground state of the system is 
formed in elongated trap at the natural value of qbf (b) 
the strength of attraction qbf is then increased. Another 
way could be performing the evaporation already under 
favorable conditions, i.e., at the presence of the appro- 
priate magnetic field or strong enough one-dimensional 
geometry. However, in this case only one, placed at the 
center of the trap double-peak soliton is formed. 

In conclusion, we have shown that bright solitons 
can be generated in a Bose-Fermi mixture as a result 
of a competition between two interparticle interactions: 
boson-boson collisions which are effectively repulsive and 
boson-fermion collisions which are attractive. Assuming 
the strength of attraction is large enough both kinds of 
atoms start to mediate in the other species interaction 
introducing the system into a new regime where locally 
Bose and Fermi gases become gases of effectively attrac- 
tive atoms. Therefore it becomes possible to generate 
bright solitons in the system under such conditions. De- 
pending on how fast the change of the attraction strength 
is performed the system responses forming the train of 
solitons (fast change) or the single soliton (adiabatic 
change). Each soliton is, in fact, the double-peak struc- 
ture with the fermionic cloud hidden within the bosonic 
one. 
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